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A mathematical model of a symmetric rotor supported by one rigid and one fluid lub- 
ricated bearing is proposed in this paper. The rotor model is represented by gener- 
alized (modal) parameters of its first bending mode. The rotational character of 
the bearing fluid force is taken into account. The model yields synchronous vibra- 
tions due to rotor unbalance as a particular solution of the equations of motion, 
rotor/bearing system natural frequencies and corresponding self-excited vibrations 
known as "oil whirl" and "oil whip." The stability analysis yields rotative speed 
threshold of stability. The model also gives the evaluation of stability of the 
rotor synchronous vibrations. In the first balance resonance speed region two more 
thresholds of stability are yielded. The width of this stability region is directly 
related to the amount of rotor unbalance. 

The results of the analysis based on this model stand with very good agreement 
with field observations of rotor dynamic behavior and the experimental results. 


1. INTRODUCTION 


Dynamic phenomena induced by interaction between the rotor and bearing or seal fluid 
motion and creating severe rotor vibrations have been recognized for over 50 years. 
High amplitude shaft vibrations which can sustain over a wide range of rotational 
speeds not only perturb the normal operation of a rotating machine, but may also 
cause serious damage to the machine and the entire plant. 

Literature related to rotor/bearing and rotor/seal phenomena is very rich [1-7]. 
Availability of computers and fast development of numerical methods are bringing 
more and more results based on analytical models of the solid/fluid interaction 
phenomena. 

There is still, however, a big gap between theory and practice. 

Practical rotating machinery instability problems in the field are being corrected 
ad hoc by applying trial and error approach with a number of measures such as 
increasing rotor radial load, modifying lubricant temperature and/or pressure, 
shortening and stiffening the shaft, or replacing bearings or seals with "more 
stable" ones. However, the only measure which corrects the causative agent is the 
recently introduced anti-swirling technique. 

Researchers and engineers do not always agree upon the physical description of the 
shaft/bearing or shaft/seal solid/fluid interaction dynamic phenomena. The complex- 
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ity of these phenomena and the long list of factors affecting them make the picture 
tremendously obscure. 

Most often bearing and/or seal fluid forces generated during rotating machine opera- 
tion are considered separately from the shaft motion, assuming only that the shaft 
rotates with constant angular velocity and has a perfect geometry. 

Practically observed phenomena indicate that this approach can be justified only for 
specific conditions, in particular, for low values of shaft rotative speed. In most 
other cases, rotor and bearings (or rotor/beari ngs/seal s ) should be considered as a 
system. The dynamic behavior of this system will then reflect the rotor/bearing 
coupled system features. 

In this paper, an attempt to build a simple rotor/bearing system model is made. The 
simplicity allows for obtaining analytical solutions and yields clear conclusions on 
how various parameters of the rotor/bearing system affect its dynamic behavior. 
The model is based on modal behavior of the rotor. The model of bearing fluid 
dynamic forces is developed through a classical approach applied to nonconservative 
mechanical systems [8]. 

The model yields results which are in perfect agreement with the observed dynamic 
phenomena of the rotor/bearing system. 


2. OBSERVED VIBRATIONAL PHENOMENA 


Lightly loaded and slightly unbalanced symmetrical rotors rotating in fluid 360° 
lubricated cylindrical bearings, exhibit the following dynamic phenomena: 

(i) When the shaft starts rotating with a slowly increasing rotative speed, 
SYNCHRONOUS (lx) lateral vibrations with minor amplitudes are observed all along 
the rotor axis. These vibrations are caused by the inertia forces of unbalance of 
the rotor. At low rotative speeds, these vibrations are stable, an impulse perturb- 
ation of the rotor causes a short time transient vibration process, and the same 
vibration pattern is reestablished (Fig. 1). 

(if) At higher rotational speeds (usually below the first balance resonance), 
the forced synchronous vibration is not the only regime of vibration. Along with lx 
vibrations, OIL WHIRL appears. Oil whirl is the rotor lateral forward precessional 
subharmonic vibration around the bearing center at a frequency close to half the 
rotative speed. In this range of the rotative speed, the rotor behaves as a rigid 
body. The amplitudes of oil whirl are usually much higher than those of synchronous 
vibrations; they are, however, limited by the bearing clearance and the fluid non- 
linear forces. With increasing rotative speed, the pattern of vibration remains 
stable. The oil whirl "half 11 frequency follows the increasing rotative speed, main- 
taining the *1/2 ratio with it. The vibration amplitudes remain nearly constant and 
usually high. At the bearing, the vibration amplitude may cover nearly all bearing 
clearance. In the considered range of rotative speed, the bearing fluid dynamic 
effects clearly dominate. The forced synchronous vibration represents a small frac- 
tion of vibration response, as the spectrum indicates (Fig. 1). 
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(iii) When the increasing rotative speed approaches the first balance resonance, 
i.e., the first natural frequency of the rotor, oil whirl suddenly becomes unstable 
and disappears, being suppressed and replaced by increasing SYNCHRONOUS vibrations. 
The forced vibrations dominate, reaching the highest amplitudes at resonant fre- 
quency corresponding to the mass/stiffness/damping properties of the rotor. The 
bearing fluid dynamic effects now yield priority to the elastic rotor mechanical 
effects. 

(iv) Above the first balance resonance speed, the synchronous forced vibrations 
decay. Again, the bearing-fluid forces come back into action. With increasing 
rotative speed, shortly after the first balance resonance, OIL WHIRL occurs again. 
The previously described pattern repeats. The width of the rotative speed region in 
which the synchronous vibration dominate depends directly on the amount of rotor 
unbalance: higher unbalance, wider is this region. 

(v) When the rotative speed approaches double the value of the rotor first 
balance resonance, the half-speed oil whirl frequency reaches the value of the first 
balance resonance — the first natural frequency of the rotor. The oil whirl pat- 
tern is replaced by OIL WHIP — a lateral forward precessional subharmonic vibration 
of the rotor. Oil whip has a constant frequency: independently of the rotative 
speed increase, the oil whip frequency remains close to the first natural frequency 
of the rotor. In this range of high rotative speed, the shaft cannot be considered 
rigid. Its flexibility, i.e., additional degrees of freedom, causes that rotor/ 
bearing\system is closely coupled. The rotor parameters (its mass and stiffness, in 
particular) become the dominant dynamic factors. The amplitude of oil -whip journal 
vibration' is limited by the bearing clearance, but the shaft vibration may become 
very high,' as the shaft vibrates at its natural frequency, i.e., in the resonant 
conditions. 

In various machines furnished with fluid lubricated bearings and/or seals, the 
above-described phenomena may take various forms, as other external factors may 
affect the system dynamic behavior. Generally, however, similar patterns can be 
expected [9]. 

Among the described dynamic phenomena, there is a clear distinction related to their 
nature: (a) rotor synchronous lateral vibration (lx) due to unbalance and (b) 
rotor fluid-related vibration. The first type is typically forced vibration. The 
rotating periodic inertia force considered "external" to rotor lateral motion causes 
the rotor response with the same frequency. The resulting motion has the form of 
the classical synchronous (lx) excited vibration. 

As there is no other external force to excite the vibrations, it is quite reasonable 
to refer to the second type of vibration, as SELF-EXCITED vibrations, occurring due 
to an internal feedback mechanism transferring the rotative energy into vibrations. 
Self-excited vibrations cannot occur in a conservative or "passive" structure, with 
no energy supply (in nonrotating systems in particular). In passive structures, the 
free vibrations following an external perturbing impulse usually have a decaying 
character, due to the stabilizing effect of damping, naturally existing in the sys- 
tem. Another situation takes place if the system is subject to a constant supply of 
energy (an "active", nonconservative structure). Well recognized are wind induced 
vibrations known as flutter. The rotating machine belongs also to this category. 
The internal energy transfer mechanism, in this case a bearing fluid involvement in 
motion, uses a part of the rotative energy to create forces, having the direction 
opposite to the damping force. The result consists of reduction, then with their 
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increasing value, a nullification of effective damping, the stabilizing factor. In 
such conditions, free vibrations do not have the decaying character anymore (effect 
of negative damping). While vibration amplitude increases, nonlinear factors become 
significant and eventually amplitudes are limited. Vibrations become periodic with 
a constant amplitude. The stable limit cycle is reached. This represents the 
practically observed case of the oil whirl and oil whip. As the oil whirl and oil 
whip occur in the system having a constant energy supply, the resulting vibration is 
referred to as self-excited. The last term is also closely related to the nonlinear 
character of the phenomena. In particular, the size of the oil whirl/whip orbits 
(limit cycle vibration amplitudes) are determined exclusively by the nonlinear 
factors in the oil bearing. 

Simple linear modelization of the rotor/bearing system can be applied as a first 
approximation in mathematical description of the observed phenomena. The linear 
model provides the spectrum of natural frequencies of the system, as these frequen- 
cies are very insensitive to nonlinear factors. The linear model also gives the 
evaluation of stability threshold, closely related to the delicate balance between 
the system's natural damping and the bearing fluid forces acting in antiphase and 
opposing damping. 

Very often oil whirl and oil whip are described as "unstable" rotor motion in a 
sense which is rather close to the terms "undesirable" or "unacceptable" rotor 
vibrations. Obviously, oil whirl/whip vibrations are highly undesirable; they dis- 
turb the machine's normal operation. The "normal operation" is related to the pure 
rotational motion of the rotor, around the proper axis and following a suitable 
angular speed. This is the only regime of motion which is required. The occurrence 
of the oil whirl/whip vibrations signifies that this pure rotational regime becomes 
UNSTABLE*, and the oil whirl /whip vibrations represent a STABLE regime. The term 
"stability" is used here in the most popular sense (following Lyapunov's defini- 
tion). The pure rotational motion (meaning zero lateral vibration) is unstable, oil 
whirl/whip lateral vibrations are stable. They exist, and any impulse perturbation 
cannot significantly modify their pattern. After a short-time transient process, 
the oil whirl /whip pattern is reestablished. 

A practical stability definition for a rotating machine is discussed in [10]. 


3. FLUID FORCES 


Derivation of fluid forces in a bearing (or a seal) is based on the consideration 
that the fluid rotation (dragged into motion by shaft rotation) plays an appreciable 
role in resulting dynamic phenomena and may, therefore, have a significant effect on 
rotor vibrations [8, 11, 12]. 

In the following presentation, an assessment of this effect is attempted in simpli- 
fied terms. It is assumed that when the shaft is rotating centered, fully developed 
fluid flow is established in the circumferential direction, that is, on the average, 
the fluid is rotating at the rate Aiu„ where uj r is the shaft rotative speed and K is 
the fluid average swirling ratio which value -is close to a half (Fig. 2). It is 

* Actually, with existing residual imbalance, the rotor pure rotational motion of 
a rotor does not exist. The only regime is forced synchronous lateral vibration 
along with the rotation. "Stability" refers then to the stability of this regime. 
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supposed that shaft lateral vibrations are small enough to make modifications of 
this pattern negligible. The flow axial component is supposed to affect values of 
the fluid forces in the x,y plane in a parametric way only, i.e. , the fluid circum- 
ferential force may proportionally increase with increasing axial flow. However, it 
is assumed that there is no feedback, i.e., the fluid axial motion is uncoupled from 
the circumferential motion (and not investigated in this paper). It is also assumed 
that the axial flow does not modify the average swirling ratio \. 

The vital assumption is that the fluid force which results from averaging the cir- 
cumferential flow is rotating with angular velocity Aw R . In rotating reference 
coordinates x , y (Fig. 2), the average fluid flow is purely rectilinear and the 
fluid force is r as follows: 

F = [K 0 +«l>i(|z r i)] z r + [D+* 2 (|z r l)] z r + M f z f (1) 

where z f = x f + jy f j = -fl U r l = Vx*+y* 

In the equation (1) K 0 , D, and are fluid stiffness, damping and inertia coeffic- 
ients respectively. <1^ and t|} 2 are nonlinear functions of the radial displacement 
| z | . It is assumed that these functions have analytical character (or at least are 
continuous, with continuous first two derivatives). Later on, as an example, the 
following nonlinear functions will be analyzed (the first symmetric term of the 
Taylor series for any nonlinear analytical functionV: 

4»i> 2 = Bi, 2 lz r I 2 , (2) 

wher,e Bj. and B 2 are positive constants. 

The average fluid force (1) has, therefore, nonlinear character. Stiffness and 
damping components of the fluid force increase with increasing journal eccentricity. 

In fixed reference coordinates x,y (Fig. 2), the fluid force will have the following 
form: 

F = [Ko-hMIzI)] z + [D-Hjj 2 (|zf)] (z-jAuj r z) + M f (2-2jAiu R z-\ 2 u>|z) (3) 

jAui-t 

where z = x+jy , |z| - V* 2 + y 2 . The equation z = z e represents transformation 
of coordinates. r 

The fluid force (3) can be presented in a classical "bearing coefficient" format: 
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The last two matrices contain nonlinear components of the fluid force. As it is 
easily noticed, the fluid force is supposed to have a symmetric character: the di- 
agonal terms are identical, the off-diagonal terms are skew symmetric. More impor- 
tant, however, is the fact that the off-diagonal terms are generated as the result 
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of rotational character of the fluid force : the tangential (or "cross") damping is 
the result of Coriolis inertia force, the tangential (or "cross") stiffness is 
generated by the relative velocity and radial damping. In addition, the radial 
stiffness Kq appearing at the main diagonal of the stiffness matrix is now modified 
by the centripetal fluid inertia force which appears with the negative sign. 

During experimental testing by applying perturbation method [13], the character of 
the fluid force expressed by the equation (3) was fully confirmed. For relatively 
large clearance- to* radius ratio, the fluid inertia force become significant, and 
modifies "damping" and "stiffness" matrices considerably. The resultive radial 
stiffness can very easily reach negative values. 

Another important conclusion relates to the "cross stiffness" coefficient, the most 
important component affecting rotor stability. This term is directly generated by 
the linear radial damping D, as the result of rotating character of the damping 
force. This term is proportional to the rotative speed iu R , i.e. , its significance 
increases with rotative speed. An immediate conclusion is* that an increase of the 
bearing damping D will not help to prevent rotor instability as the "cross stiff- 
ness" increases proportionally to D. The only help in decreasing the "cross stiff- 
ness" term is a modification of the average swirling ratio k. This feature is now 
being widely used in "anti -swirling" devices [14, 15]. 


4. MATHEMATICAL MODEL OF A SYMMETRIC FLEXIBLE SHAFT ROTATING IN ONE RIGID 
ANTIFRICTION BEARING AND ONE OIL 360° LUBRICATED CYLINDRICAL BEARING 

Oil whirl and oil whip phenomena are characterized by low frequency, that is why it 
seems reasonable to limit the rotor model to its lowest bending mode. It is obvious 
that the considerations presented below can be applied to more complex rotor model 
(including more modes, gyroscopic effect, internal damping, etc.) as well. 

The mathematical model representing balance of forces in the symmetric rotor shown 
in Fig. 3 is as follows: 

MiZt + D s Zi + (K 1 +K 2 )z 1 
M 2 2 2 + M f (z 2 -2jAw R z 2 -X 2 w|z 2 ) + [D-h|j 2 (|z 2 | )](z 2 -jAui R z 2 ) + [Ko+i|>i(|z 2 |)]z 2 + (5) 

+ K 3 z 2 + K 2 (z 2 -Zi) = 0 , z x = x t + jyi , z 2 = x 2 + jy 2 

where M t , M 2 are rotor generalized (modal) masses, D is external generalized 
(viscous) damping coefficient, Kj, K 2 , and K 3 are shaft generalized (modal) stiff- 
ness coefficients (K 3 may also include an external spring stiffness), m and r are 
mass and radius of modal unbalance respectively. The bearing fluid dynamic force is 
introduced to the equation (5) in the form (3). 

The equation (4) presents the classical relationship for an unbalanced symmetric 
rotor at its first bending mode. Any classical method of modal reduction can be 
applied to obtain the rotor generalized (modal) coefficients. The equation (5) 
describes balance of forces, including nonlinear fluid force, acting at the journal, 
Radial forces (such as gravity) are balanced by external spring forces and are not 
included in the model (4), (5). 


- K 2 z 2 = mnu|e 


ju) R t 


(4) 
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5. SYNCHRONOUS SOLUTION: lx ROTOR FORCED VIBRATIONS DUE TO UNBALANCE 

In this section, the particular solution of equations (4) and (5) will be discussed. 
This solution describes rotor synchronous vibrations due to unbalance force and has 
following form! > „ . (6) 


where the amplitudes A x , A 2 and phase angles at, a 2 can be found from the following 
algebraic equations resulting from (4), (5) and (6): 

(Kt+Ka+Dgj^-Mw^Aje^ 1 - K 2 A 2 e Ja2 = mm 2 

{ Ko+K 2 +K 3 +M*x ( A 2 )+j [D+M i 2 ( A 2 ) 3ui R (l-A) - M 2 u»j-M f u ) 2(l-X)2}A 2 ej“ 2 -K 2 A 1 e Jai = 0 

Note that the fluid Inertia and radial damping carry coefflcents 1 - x, which value Is 
usually slightly higher than 1/2. This signifies that the effective damping Is two times 
and the fluid Inertia effect on synchronous vibrations four times, smaller than their 
actual values. 


The equations (7) can be solved for any given nonlinear functions il^, and ij> 2 , as 
their arguments now become equal to A 2 , the journal amplitude of synchronous vibra- 
tions. When the nonlinear forces are neglected («j>i=»|<2 = 0) the equations (7) yield 
simple expressions for the amplitudes and phase angles: 
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"1 = arctan <-0 s V h '> ) + arctan l^AMaeos 'Sn ■ “2 = arctan (-h,/h,) (9) 

or, further R 

mrtUn 

A t = V(h 2 +Kl) z +h§ , ai - arctan (-D_u)n/tu) + arctan [K2h a /(h$+h<H< 2 h 2 )] 

Vh|+D|w| s K (10) 

where h 2 = h 4 h 5 - DD s u)|(l-\)-Kl , h 3 = [h s D s +Dh 4 (l-\)]u) R (11) 

h 4 = Kt+Ka-MiuJ , h 5 = K 0 +K 2 +K 3 -u)|[M 2 +M f (l-\) 2 ] (12) 

When the nonlinear forces are given in the form (2), then the amplitude A 2 can be 
found as the solution of the following polynomial equation: 


Af(h§+h?) + 2A|(h 2 h 6 +h 3 h 7 ) + Al(h§+h§) - (K 2 mnuj) 2 = 0 (13) 

where h 6 = Bih 4 - u»2B 2 (l-\) , h 7 = w R [B 2 (l-\)h 4 +B 1 D s ] (14) 

The phase a 2 in this case is slightly modified by the nonlinear force: 

Aiw R [D B 1 +(l-\)h 4 B 2 ]+h 3 

" a = arctan [ ^F R 7 1->7i3 5 B ; ,-h 4 B 1 ]-h ; , J (15) 

The amplitude A t , and phase a x are given by equations (8) and (9). 
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The equation (13) can be solved graphically (Fig. 4). The simple graphical method 
reveals important qualitative features of the nonlinear system. It is possible that 
three solutions for A 2 are yielded (Fig. 4a). This corresponds to. the well known 
nonlinear case [16], More probable, however, is that only one solution, exists. The 
nonlinear force causes a reduction of the vibration amplitude A 2 . Figure 5 gives 
some numerical examples of Bode plots for the rotor synchronous response. An 
insensitivity of the response to fluid inertia and journal generalized (modal) mass 
is noted. 


6, STABILITY OF THE PURE ROTATIONAL MOTION OF THE SHAFT 


To investigate the stability of the pure rotational motion of the shaft, assume that 
the unbalance force is equal to zero. For the stability analysis, the nonlinear 
terms in the equation (5) are neglected, as they are small of the second order. 
Stability of the pure rotational motion of the shaft means stability of the zero 
solution of the equations (4) and (5) for mr=0, perturbed by small lateral deflec- 
tions. The eigenvalue problem with the solution for linearized equations (4), (5) 
of the form . . 

Zi Eie , z 2 ~ E 2 e (16) 

leads to the corresponding characteristic equation: 

[M 2 s 2 + M f (s-jAu) R ) 2 + D(s-jAiu R ) + K 0 +K 3 ] ( Kj+Ka+DgS+Mi s 2 ) + M^+DgS+MiS 2 ) = 0 (17) 

where s is the system eigenvalue, and E lf E 2 are constants of integration. The 
equation (17) can easily be solved numerically. An example is presented in Fi gure 
6. The eigenvalues are gi ven as functions of rotative speed ui R . 

After computi ng several examples, it has been noticed that the modal mass M 2 and the 
f 1 ui d i nerti a M~, when not exceeding certai n "critical" values, have very little 
effect on the eigenvalues. It has also been noticed that one of the eigenvalues has 
always the imaginary part proportional to the rotative speed u» R , when the latter is 
low, an d that it t ends to the constant value corresponding to tne rotor natural fre- 
quency V(Ki + K 2 )/M t when uj r increases. The corresponding imaginary part of this 
eigenvalue crosses zero dt a specific rotative speed. That is, after this specific 
speed, the pure rotational motion of the rotor becomes unstable. 

These observations led to an important conclusion about the character of the system 
eigenvalues, and to approximate eigenvalue formulas. 

The approximate values of three eigenvalues, 1.e. f three solutions of equation (17) are 
as follows (obtained from entenslve parametric analysis of numerical solutions to 
eq. (17)): 


+K „ MKi-MiA 2 * 2 ) 

IT* " D(K 1 +K 2 -M 1 \ 2 uj 2 ) + jAu) R 


(18) 
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s 2 ,3 * ± jVRpjRl = ± jVRi + >/RFRI)/V2 


(19) 


where 


R.i * - K|(Ko+K 3 )/R 3 , Ra * -KiD(V^^-Au JR )/R 3 

R 3 = {[K 0 +K 3 -M 2 (K 1 +K 2 )/M 1 - M f (V(K 1 +K 2 )/M 1 -\u) R )2]2 + D2(V^^-A 10 r )2}M 1 


The formulas (18), (19) give approximate eigenvalues of the rotor/bearing model (4), 
(5). The imaginary parts of the eigenvalues (18) and (19) represent natural fre- 
quencies. The first natural frequency, Im (s*) is close to Xm R . The seco nd natural 
frequency, Im (s 2 , 3 ) is close to the shaft natural frequency ±" V(Ki + K 2 )/Mi. 

Stronger coupling (high K 2 ) causes more significant divergence from these 
"uncoupled" natural frequencies of the system. 


The real part of the eigenvalue (18) predicts the threshold of stability: 


“R i 
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( 20 ) 


the rotor pure rotative motion is stable. For u»„ > iu R 'the pure rotational motion 
becomes unstable. The first term under the radical *20), mainly Ki/Mj, is defin- 
itely dominant. The second term contains two stiffnesses in sequence, K 2 and 
(Ko+Ks-I^Ki/Mi). Usually the stiffnesses K 0 and K 3 are small, mass M 2 is also 
small. Connection in sequence with K 2 (independently of value K 2 ) makes the total 
smaller than K 0 +K 3 -M 2 K 1 /M 1 . It is reasonable, therefore, to reduce the expression 
for the rotor/bearing system stabi 1 i ty threshold to 


U) 


(ST) 

R 



( 21 ) 


For the unstable conditions the vibration amplitude increases exponentially in time 
and eventually bearing nonlinear forces become significant causing final limitation 
and stabilization of the vibration amplitude. 


The experimental results entirely confirm this result. An increase of the stiffness 
Ki (disk mounted on the shaft moved toward the rigid antifriction bearing) causes a 
significant increase of the system threshold of stability [17]. 


7. STABILITY OF THE SYNCHRONOUS SOLUTION 


The stability of the synchronous vibrations (6) will be analyzed by applying the 
classical perturbation, method [16]. Introducing the variational real variables 
w x (t), w 2 (t) (amplitude perturbation) according to the relation 

Zi = [A 1 +w 1 (t)]e j(u, R t+ofl) , z 2 = [A 2 +w 2 (t)]e J ’ (u, R t+a2) (22) 

the linearized variational equations are obtained by substituting Into equations (4) and 
(5). 
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Mi ( w x +2w j.u» R J w x ) + (K 1 +K 2 )w 1 -K 2 w 2 e j(a2 " ai) = 0 
M 2 (w 2 +2w 2 uu R j-iJu|w 2 ) + M f [w 2 +2ju) R w 2 (l-\)-u)|w 2 (l-\) 2 ] + 

+ [D+4» 2 (A 2 )] [w 2 +jw R w 2 (l-\)] + w 2 A 2 juj r (1-\) + 

+ CKo-N»i(A 2 >] w 2 + w 2 A 2 + (K 2 +K 3 ) w 2 - K 2 w ie j(ai " a2) = 0 

The solution of the variational equations (23) has the following form 

w, = E x e (st ‘ j “ , R t ‘ J '“ l) , w 2 = E 2 e (st '*"R t " J “ 2) (24) 

where s is the eigenvalue and E*, E 2 are constants of integration. By introducing 
(24) into (23), the characteristic equation is obtained: 

[M 2 s 2 +M f (s-jAw R ) 2 + (D+t}> 2 )(s-jAu) R ) + »J» 2 A 2 ju) R (l-\) + 

+ K 0 -H|ii-H|»iA 2 +K3](K 1 +K 2 +O s s+M 1 s 2 ) + K 2 (Ki+D s s+MiS 2 ) = 0 

where tji! and <J> 2 are functions of A 2 and ij>i, 2 = • 

Stability of the synchronous solution (6) is assured if all eigenvalues of the equa- 
tions (23) have non-positive real parts. 


The characteristic equation (25) di f fers from the characteristic equation (17) by 
the nonlinear terms which cause an apparent i ncrease of the beari ng radi al damping, 
"cross stiffness" effect, and radi al stiffness. The threshold of stability cal cu- 
lated the same way as previously, has now the following form (M 2 is neglected): 

CST> — 1 — a, * mjsa p (26) 

WRTCi + x( a 2 )] 


UJn — 


where x(A 2 ) = -<|» 2 A 2 (l/\-l)/(D+t|» 2 ). 


(27) 


With increasing eccentricity A 2 an apparent decrease of the average swirling ratio k 
is noticed, as x(A 2 ) has a negative sign. 


More interesting is, however, the modification of the instability threshold due to 
increase of stiffness from (K 0 +K 3 ) in (20) to (K 0 +K 3 -h|j 1 -hjjiA 2 ) in (26). The last 
expression includes the amplitude of synchronous vibration A 2 which is a function of 
rotative speed tu p . Following (20) and (26), the criterion of stability can be 
wri tten in the for 


fi » i v wwii I i i vi iw i v t • 

< K 7 .yK, + 


K 2 \|i 


(28) 


and solved graphically (Fig. 7). The graphical solution is very effective in show- 
ing qualitative features of the rotor dynamic behavior. Having the amplitude/ 
rotative speed relationship from equation (8) or (13) (lower graph of Fig. 7) and 
having the plot of the radical (28) versus amplitude A 2 which includes the nonlinear 
function of A 2 (left-hand side graph), the graph of the radical from the right- 
hand side of the inequality (28) versus rotative speed uj r can be built up as the 
small arrows indicate (right side plot). The left-hand side of the inequality (28) 
versus m R is a straight line Au) R -/M t slightly modified by the factor x(A 2 ). The 
regions ef stability are founcT at intersections of the corresponding plots on the 
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upper right-hand side graph. It clearly shows that the synchronous solution is 
stable for rotative speeds below the first threshold of stability, which differs very 
little from the stability threshold (21), and that they are stable in the resonant 
speed region. 

It is easily noticed that higher unbalance will i ncrease sta bi 1 i ty region for the 
synchronous vibrations around the resonant speed V(Kx + K 2 )/M 1 (Fig. 8). 

The series of experiments with a balanced and unbalanced rotor confirm this analyti- 
cal prediction; higher unbalance causes wider interval of rotative speeds where 
synchronous vibrations are stable (Figs. 9 through 11). 


8. ROTOR SELF-EXCITED VIBRATIONS: OIL WHIRL AND OIL WHIP 


For the case of the absence of imbalance, the equations (4), (5) have a particular 
solution describing self-excited vibrations (oil whirl and oil whip). 

Zi = Gl e J * (u)t+P) , z 2 = G 2 e JU,t (29) 

where G t , G 2 are corresponding constant amplitudes, p is the relative phase angle, 
and u» is the self-excited vibration frequency. 

By introducing (29) into (4) and (5), the algebraic equations for calculating G lt 
G 2 , p, and w are obtained: 

{-M 2 u)2-M f (u»-Au) R )2 + j[D+iJj 2 (G 2 )](ur\u)p) + Ko+<h(G 2 ) + K 3 ](K 1 +K 2 -M 1 u)2+D e juj) + 

K S (30) 

+ K 2 (Ki+D jui-MiU) 2 ) = 0 


Gi = 


K?G? 


D s u» 


4 — 4 , P = arctan — t-tt — — 

V(Ki+K 2 -M lU )ii) + D*u»* 


(31) 


The real and imaginary parts of the equation (30) provide relationships for obtain- 
ing G 2 and ui, for each case of the nonlinear functions $1(62) and 0 2 (G 2 ): 


(K 1 +K 2 -M 1 uj2)(K 1 -M 1 u) 2) + D|iu 2 


*!>i(G 2 ) - M 2 u)2+M f (u)-Au) R )2-K 0 -K 3 -K 2 (K i+K 2 -M 1 u)2)2 + ^ 


= <j>l(u>) 


^2(^2) - -D - 


K1D s uj 

(u»-Am R )[(K 1 +K 2 -M 1 uj2)2 + D|u»2] 


= 02 ( u) ) 


(32) 

(33) 


The frequency w is yielded from the equation 

4>2{0 i - 1 l 0 i(w)]} = 0 2 (uO (34) 

where 0! -1 is the inverse function (32); G 2 = 4 ji _ 1 [0 i(uj)]. 

For example, in case of the "parabolic" nonlinear functions (2) and small damping D g 
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the frequency equation (34) has the following form: 
(K 1 +K 2 -M 1 uj2)(u J -Aa) R ){(K 1 +K 2 -MiU)2)[|^(Ko+K3-M 2 u)2-M f (u)-Au) R )2)-D] + 

+ §^K 2 (Ki-MiU)2 )} * D s u)Kl * 0 

It yields immediately three important real solutions for the frequency: 

UJj * Auu r 


u ) 2 , 3 * ± V(Ki+K 2 )/Mi. + D s Bi/[M 1 B2(l-Au» R /V(Ki+R 2 )7MrrT] * ±V(Ki+K 2 )/Mi (37) 

which are close to the linear rotor/bearing system "whirl" and "whip" natural fre- 
quencies at the threshold of stability (21). The remaining solutions can be also 
easily obtained from equation (35). For the small shaft damping D the similar 
consideration can be extended to. any case of nonlinear functions <1 »i, s i)j 2 : in (35) 

the ratio B 2 /Bi should be replaced by the ratio «|» 2 (G 2 )/t|» 1 (G 2 ). Three frequency 
solutions (36), (37) exist, therefore, independently of nonlinearities and independ- 
ently of journal mass and fluid inertia. 

The corresponding amplitude G 2 of the self-excited vibrations can be calculated from 
(32). For example, the "parabolic" nonlinearity (2) yields: 


( 35 ) 


(36) 


/ Kt-MiXZujS 

Whirl amplitude: G 2 | » ^ Ax^-Ko-Ka-K, K.+Kg-M.^ 


(38) 


Whip amplitude: G 2 | 


uW(Ki+K 2 )/M 


* 3 - VM2(Ki+K 2 )/M 1 +M f (V(K 1 +k 2 )/R7-Au )R )2 - K 0 -K 2 -K 3 

(39) 


The expression under the radical (38) should be positive, therefore, the self- 
excited vibrations with the whirl amplitude (38) exist only for the limited range of 
the rotative speed (M 2 neglected): 


l "fe + ferfe+K 2 +K 3 ) ~ W R ~ V(l<i + K 2 )/M 1 A ( 40 ) 

The left-hand side of the inequality (40) represents the threshold of stability 
( 20 ), the right-hand side term separates the whirl from the whip. 


The similar reason causes that the self-excited vibrations with the whip amplitude 
(39) exist for the following rotative speed range: 


u> R > ± i V C Ko+K 2 +K 3 - ( K x +K 2 )M 2 /Mi J/M f 


(41) 


From equations (31) and (32) the corresponding whirl and whip amplitudes and phases 
of the rotor can be obtained. For example for whip frequency (37), in case of 
"parabolic" nonlinearities ( 2 ) they are as follows: 
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G x = G 2 


( 42 ) 




0 <; V(Ki+K 2 )/Mi Vb!/[B 2 (V(Kx+K 2 )/Mx-Au) r )]2 + 1 

u)=V(Kx+K 2 )/M s * 


P = arctan [B 2 (V(K 1 +K 2 )/M 1 "Ajuu r )/B 1 ] (43) 

During whip, the shaft vibrates at its resonant conditions, the whip amplitude is 
controlled mainly by the shaft damping D g and the nonlinear fluid forces. 

The phase between the rotor and journal at whip vibrations is controlled by non- 
linear forces in the bearing. It was experimentally verified that during both oil 
whirl and oil whip, the rotor and journal vibrate nearly in phase (Fig. 12). This 
may serve as an indication concerning the values of nonlinear forces t|>x and ij» 2 . 


9. DISCUSSION OF RESULTS 


In this paper, a simple mathematical model of a symmetric shaft rotating in one 
rigid and one 360° lubricated bearing is proposed. The model yields results which 
stand in very good agreement with the experimentally observed rotor dynamic phenom- 
ena concerning stability, thresholds, and the self-excited vibrations known as oil 
whirl and oil whip. 

The rotor in the model is represented by the generalized (modal) parameters of its 
first bending mode. These parameters can be analytically obtained by applying any 
classical method of modal reduction. The fluid force acting at the journal is 
represented in the model by the nonlinear expression based on rotating character of 
this force. The latter approach yields the bearing coefficients in an "inter- 
related 11 form: the tangential ("cross") stiffness is a product of the relative vel- 
ocity and radial damping, the tangential damping results from Coriolis acceleration, 
the radial stiffness is reduced due to centripetal acceleration and fluid inertia. 

The fluid nonlinear force is introduced in a very general form, where the non- 
linearity is related to the journal radial displacement (eccentricity in the journal 
bearing). The nonlinear character of the fluid stiffness force and especially fluid 
damping force is well known. This nonlinearity has a "hard" character: the forces 
increase with eccentricity. Actually they grow to infinity when the journal touches 
the bearing wall. Since, in this paper, low eccentricity ratios were considered 
only (to justify the symmetric character of the rotor model), an example of the 
fluid nonlinear forces was discussed: the first terms of the Fourier expansion of 
any nonlinear function representing the fluid force.* 

The results obtained from the analysis of the rotor model and concerning the rotor, 
synchronous vibrations due to imbalance, and self-excited vibrations known as "oil 
whirl" and "oil whip" very well reflect the observed rotor dynamical behavior. 

The classical eigenvalue problem yields three important eigenvalues of the rotor/ 
bearing system. The first eigenvalue has the imaginary part (natural frequency) 
close to and corresponds to oil whirl frequency. The real part of this eigen- 
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value predicts the threshold of stability (the rotative speed at which the pure 
rotational motion or lx synchronous vibrations of an unbalanced rotor become 
unstable). This threshold of stability is determined by the rotoE-p m ass an d its 
partial stiffness, as well as the oil average swirling ratio (u) R '- J ~ 
expression (21)). The latter result was noticed in the fifties oy Poritsky (18] and 
discussed by Boeker and Sternlicht [19]. Since then, this seems to have been 
forgotten for 30 years. 

The second and third eigenvalue of the rotor/beari ng system have the imaginary part 
close to the rotor natural frequency (±V(K 1 +k 2 )/M 1 ) of the first bending mode. The 
latter is referred to as "whip frequency". 

The model yields the self-excited vibrations (known as oil whi rl and oil whip) as a 
particular solution. The frequencies of oil whi rl and oil whip are very close to 
the natural frequencies of the 1 i near system. The model permits evaluation of the 
amplitudes and relative, journal /disk phase angles of these self- excited vibrations. 

The most important and new (to the author 1 s knowledge) result presented in this 
paper concerns the analytical evaluation of the stability of the synchronous vibra- 
tions of the rotor and relationship between the width of the stability region and 
the amount of the rotor imbalance. This result is yielded from the classical 
investigation of stability through equati ons in variations. The beari ng non! inear 
force, increasing with rotor radi al def 1 ecti on causes an increase of the journal 
supporting radial stiffness force which is directly related to the actual synchro- 
nous vibration amplitude. The i ncreased radial stiffness affects the threshold of 
stability. The synchronous vibration amplitude varies with the rotative speed, and 
in the first balance resonance region of speeds the effect of stable synchronous 
vibrations is observed. Two additional thresholds of stability are noted. The 
width of this region depends di rectly on the amount of imbalance in the rotati ng 
system. A simple graphical method (Figs. 7 and 8) explains qualitatively this 
feature. The expressions (26) and (28) give the analytical relationships for these 
additional thresholds of stability. 

In the rotor/beari ng system analysis, as well as in experimental testing, rather 
small effect of the beari ng fluid inertia has been noted. The f 1 ui d i nerti a has a 
negligible influence on the synchronous vibrations, thresholds of stability, natural 
f requenci es , and frequenci es of the self-excited vibrations. The journal general- 
ized (modal) mass is also usually relatively small. These aspects led to the con- 
clusion regardi ng further simpl ifications of the rotor/beari ng model to the one 
and a half degree of freedom system. The fluid radial damping force which is pro- 
port Tonal to the relative journal velocity is definitely domi nant fluid force in 
the bearing, and it determi nes the order of the equation (5), as the forces of 
inertia M^. and M 2 are negligible. 

There is sti 1 1 no clear picture of what affects the fluid average swi rl i ng ratio X. 
It has been assumed in this paper that X is constant. However, during experimental 
testi ng some variations of value of X have been noticed. This aspect is bei ng 
studied and requi res further investigation. 


* Fluid force is often modeled proportional to (1-e 2 ) n , where e is the eccentricity 
ratio (journal radial displacement |z 2 | to radial clearance), n is a number (usually 
1/2, 1, 3/2, 2, 5/2, or 3 ). In the first approximation (1-e 2 ) ~ 1 + 2ne 2 , which 
in turn, is proportional to the example forces (2) considered in this paper. 
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SYMBOLS 


A x , A2 

Bx , B 2 
D, Ko, M f 

D s 

Ei, E, 

F 

Gi , G 2 

j = /l 


Amplitudes of synchronous vibrations of the rotor and journal 
respectively 

Coefficients of bearing nonlinear "parabolic" force (2) 

Bearing fluid radial damping, stiffness, and inertia coefficients 
Shaft external damping coefficient 
Constants of integration 
Bearing fluid force 

Amplitudes of the rotor and journal self-excited vibrations 


Kt, K 2 , K 3 Rotor modal stiffness coefficients 

m, r Mass and radius of imbalance respectively 

Mj., M 2 Rotor and journal modal masses respectively 

s, Si, s 2 , s 3 Eigenvalues 

wx(t), w 2 (t) Variational variables for the synchronous vibration/stabil ity 
investigation 

z=x+jy, z 2 =x 2 +jy 2 Journal radial displacement x-horizontal , y-vertical (fixed 
coordinates) 

Zi = Xx+jyx Rotor radial displacement (xi-horizontal , yi- vertical) 

z_ = x_+jy_ Journal radial displacement (rotating coordinates) 

Rotor and journal phase of the synchronous vibrations 
Phase of rotor/journal self-excited vibrations 
Fluid average swirling ratio 

Bearing fluid nonlinear functions of journal radial displacement 
Derivatives of 4>i. 1° respect to radial displacement 

correspondingly 

w Circular frequency of the rotor self-excited vibration's 

Rotative speed 


P 

\ 

«l»i. 


= * r +jy r 

. ct 2 


^2 


UJ 


R 

(ST) 


Threshold of stability (rotative speed). 
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mounter uventsmmm i ioogi 

Figure 1. - Cascade spectrum of rotor vibrational response measured at oil bearing. Spectrum 
indicates regions of synchronous (lx) vibrations due to unbalance, oil whirl («1/2X), and oil 
whip with constant frequency close to first balance resonance (slightly lower), (a) Orbit display 
of synchronous vibrations, (b) oil whirl orbit, (c) oil whip orbit (K<J> - keyphasor pulse 
indicating relationship with rotative speed frequency). 



Figure 2. - Bearing mode: Fluid force rotating with angular velocity Xwr. 
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z i = *1 + jy t 


z 2 - *2 + jy 2 


Figure 3. - Model of symmetric rotor in one antifriction and one oil bearing. 


A| (a| b, +2A| b 2 + bjlJ 





LINEAR 

| | SOLUTION (8) 

1 I / 


(a) Case 2b^ < — y3b^b^: three amplitudes yielded, 
(b) Case 2b ? > one amplitude A,, exists. 


2 2 

Figure 4. - Graphical solution of equation (13) (notation: b, = h, + h,, b, = h_h, + 

p a I O / £ 4 0 

b„ = h„ + hf). 



0 1 2 3 4 

ROTATIVE SPEED CO* [RPMX 1000] 



ROTATIVE SPEED a)* IRPM X 1000] 


Figure 5. - Phase and amplitude of rotor and journal synchronous vibrations (6) versus rotative 
speed. Note a significant insensitivity to fluid inertia and journal generalized mass 


173 








Figure 8. - Hodification of synchronous vibration stability regions by unbalance same as in 
fig. 7). 



Figure 9. - Cascade spectrum of rotor with residual unbalance. Vertical response measured at 
journal. Synchronous vibrations unstable for rotative speed higher then ~2100 rpm. 
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Figure 10. - Cascade spectrum of rotor with 1 g x 1.2 in. of imbalance. Vertical response of 
journal. Synchronous vibrations stable for rotative speed lower than ~2100 rpm and in resonance 
region: from ~4200 to ~490Q rpm. 



Figure 11. - Cascade spectrum of rotor with 1.2 g x 1.2 in. of imbalance. Vertical response of 
journal. Synchronous vibrations stable for rotative speed lower than ~2100 rpm and in wider 
resonance region: from ~410Q to ~5200 rpm. 
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Figure 12. - Vertical transducer signals of journal and rotor vibrations during oil whip. 
Oscilloscope display in XV mode indicates nearly in-phase vibrations. 
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